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Abstract 


An integral equation approach, is developed to determine the scattering 
and absorption of electromagnetic radiation by thin walled cylinders of arbitrary 
cross-section and refractive index. Based on this method extensive numerical 
data are presented at wavelengths in the infrared for hollow hexagonal cross 
section cylinders which simulate columnar sheath ice crystals. 


1. latroductioE 

This is the Final Repoii; on NASA Grant NSG 5044 and describes the work 
carried out in computing the scatteriug and absorption of electromagnetic radiation 
by columnar sheath ice crystals. In this first section we begin by citing some of 
the important geophysical problems which require such data, and note some of the 
consequences of the lack of sufficiently accurate values for the scattering and 
absorption properties of single ice crystals. After surveying the work done in 
recent years to remedy this deficiency, we then summarise our own contributions 
to the subject. The detailed results are given in the remaining sections of the 
report. 

The reflection, transmission and absorption of visible and infrared radiation 
by clouds and by polluted atmospheres is of considerable importance in many prac- 
tical areas. Cirrus clouds are composed of ice crystals. They are found over the 
entire globe and their infrared and optical scattering properties have a profound 
effect on the atmospheric heat balance. The clouds s catter and absorb the primarily 
infrared radiation resulting from the earth and the lower regions of the atmosphere as 
well as the mainly shorter wavelength solar radiation incident from above. The differ- 
ence then plays a major part in the atmospheric heat balance which governs the global 
location of energy sources and sinks and, hence, the atmospheric circulation patterns 
(Cox, 1971; AFWG, 1972). Similarly, but on a much smaller scale, the cirrus clouds 
created artificially by jet contrails have been observed to markedly affect local 
weather (Appleman, 1966; Reinking, 1968); and a knowledge of the scattering by 
clouds is also needed in using satellite measurements of the IR emission of water ■ 

vapor to estimate the relative humidity of regions above the clouds. Finally, we > 

remark that the use of LIDAR as an atmospheric problmg tool depends on the dif- i 

ference in the reradiation of the aerosols and the much smaller baclcground mole- 
cules (Grams, 1975). 
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Techniques for calculating the transfer of electromagnetic radiation through 
clouds of particles have been summarized by Plass (1973), and one of their 

most basic ii^redients is a knowledge of the scattering and absorption properties 
of the individual particles. This may involve either single scattering or multiple 
scattering in the case of optically thick clouds. Since each scattering event can 
affect the polarisation by producing’ an electric field having a component orthogonal 
to the incident vector as well as parallel to it, and this field is in turn incident on 
another particle, an accurate treatment of the transfer problem will involve the com- 
plete scattering and absorption matrices for a single particle. 

In an atmospheric cloud, the particles are either water droplets or 
ice crystals . The water drops are close to spherical and it is not unreasonable to 
model them as spheres. This enables Mie theory to be applied and the results 
obtained are reasonably accurate. Ice crystals, however, are another story. The 
shapes, sizes, concentrations and fall patterns which have been observed in clouds 
have been discussed by Mossop and Ono;, 1969; Oao, 1969; Aufm Kampe and 
Weickman, 1957; andHe;^sfield and Knollenberg, 1972. Both plate crystals, i.e. 
cylinders of length much less than their diameter, and columnar crystals, i. e. 
long thin cylinders, hollow as v'ell as solid, are commonly found, and for shapes 
as varied as this, a sphere cannot provide an accurate simulation of the scattering 
behavior. Nevertheless, for lack of a more accurate method for calculating the 
scattering properties, the Mie theory has been widely used even for ice crystals, 
thereby introducing unlmown and possibly large errors in the values for the radia- 
tion transfer, which are the end products of extensive and expensive computations 
(Kattewar and Plass, 1972). 

The importance of using the proper scattering matrix when the particles are 
irregular is clear from the data presented by Holland and Gagne ( 1970) . They mea- 
sured the elements of the scattering matrix for clouds of irregularly shaped, random- 
ly oriented silicon flakes. The results were quite different from the matrix elements 
predicted by Mie theory, particularly for back and forward scatter, and the depolari- 
zation was also poorly predicted by the theory. 
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The last few years have seen, several attempts to calculate the scattering 
of more realistically-shaped crystals, and it is appropriate to mention here the work 
of Jacobowitz (1971) and Liou (1972 a and bj 1973) which has been directed at the 
scattering properties of columnar ice crystals. Jacobowitz's data were obtained for 
infinitely long crystals hexagonal in cross section using ray tracing. Ail end effects 
were necessarily omitted, including the 45° deviation of the rays passing through the 
end feces which contributes to the large halo observed about ice clouds (see Minnaert, 
1954, section 104). The method also excludes all diffraction effects produced, for 
example, by the six longitudinal edges of the cylinder, as well as polarization effects, 
and the calculations were limited to cylinders not less than 40p in diameter (for a 
wavelength of 0.55 i^) with the apparent objective of assuring the reasonable validity of 
geometrical optics. Finally, no account was taken of internal absorption by the ice 
in spite of the fact that the appreciable imaginary part of the refractive index at some 
infrared wavelengths suggests that absorption may not be negligible. 

Lion's analyses are based on the assumption that the ice crystal can be model- 
led by an infinitely long, homogeneous dielectric cylinder of circular cros^ section. 
For this simplified geometry there is a mathematically exact expression for the 
scattered field in the form of a series of Bessel and Hankel functions analogous to 
the Mie series for a sphere. It is therefore possible to compute the scattering matrix 
precisely, with all polarization information present, and with internal absorption 
taken into account. Nevertheless, end effects are omitted by virtue of the model 
chosen, and the assumption of a circular cylinder necessarily suppresses those 
features of the scattering which are peculiar to the hexagonal cross section of an 
actual ice crystal. 

The retention of the hexagonal geometry is one of the key features of the work 
carried out under the present Grant. Based on an integral equation approach, a 
numerical technique has been developed to compute the scattering patterns, and the 
scattering and absorption spectra for cylindrical dielectric shells of arbitrary cross 
sectional shape when illuminated by a plane wave of either principal polarization. 

The dielectric can be lossy, and by applying the method to infinitely long cylinders 


a 


hexagonal in cross sectionj scattering and absorption data have been generated appli- 
cable to hollow columnar {sheath) ice crystals in the infrared. 

The method originated from a study of the scattering properties of resistive 
sheets and membranes (Knott and Senior, 1974) in which the non-zero thickness 
sheets were simulated by infinitesimally thin sheets of appropriate electromagnetic 
properties. Accordingly, a hexagonal shell cylinder whose actual walls are composed 
of a material of (complex) dielectric constant n, having thickness t( small mmpared 
to the free space wavelength X), is replaced by a hexagonal membrane having a com- 

s 

plex relative resistivity 

_R ^ a 
^ (n^-l)2iTT 

ohms per square. It is then possible to derive an integral equation for the current 
which an incident plane wave of either principal polarizatiorf induces in the membrane, 
and the integral equation is quite amenable to solution by digital techniques. The 
formulation of the equation and its subsequent solution constitute a significant con- 
tribution to the theory and application of integral equation methods to electromag- 
netic scattering and absorption problems. The mathematical details are given in the 
Appendix along with a listing of the computer programs employed in generating the 
data in this Report. For the most part the formulas are applicable for arbitrary 
angles of incidence, but the numerical results presented here are for broadside 
incidence only. The problem is then two-dimensional, and we now turn to a presen- 
tation and discussion of the data obtained. 
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Nitraerical Eesults 


Our numerical results are given in th.e form of cross sections which are 
defined as follows . For a power density S incident on the cylinder, the bistatic 
scattering cross section is 

O<0) = 27TI/S 

where ;i is the power scattered per unit length of the cylinder per unit angle about 
the direction 0 and measured in the far field of the cylinder. The angle d is defined 
so that 0 = 0 is in the backscattering direction and 6 = v is in the forward. The 
total (or integrated) scatterir^ cross section is then 


CT 


T 


2sr 



o(0)d0 . 


The absorbed power is measured by the absorption cross section 

a . =* "T (power absorbed) 

A o 

and the extinction cross section is the sum 


ext 


= cr_ + a , 
T A 


Formulas relating these two dimensional cross sections to the currents which are 
computed are given in the Appendix, eqs. (37) through (41/, 

The computations were carried out using the refractive indices n = n^ + in. 

for ice in the infrared wavelength range given by Irvine and Pollack ( 1968; hereafter 
referred to as IP) and Sehaaf and Williams (1973; referred to as SW). Their data 
are plotted in Figs. 1 and 2 and show significant discrepancies in certain wave- 
length ranges. To obtain some idea of how’ sensitive the scattering is to the partic- 
ular refractive index chosen, computations have been made at two wavelengths using 
the values from each reference. 

For a given wavelength and perimeter of the hexagonal sheath, the 
scattering has been computed for the incident plane wave polarized with its electric 






V6ctor p 3 .rall 6 l to th© axis of th© cylindsr {E polarization) and. also with its magnotic 
vector parallel to the axis (H polarization). The directions of incidence and obseiwa- 
tion are always in a plane perpendicular to the axis, but for each polarization we have 
considered two directions of incidence corresponding to 'edge-on' and 'face-on' as 
regards the hexagon, viz. 



edge- on fe.ce-on 


Most of the calculations have been for a hexagonal cylinder 3 on a side 
with a wall thickness t = 0. 1 /urn. Only these data are presented here though we have 
carried out some exploratory calculations for other parameter values. Table 1 lists 
the wavelength, the corresponding refractive index and its source, the appropriate 
resistivily value and the Figure numbers where the computer-generated plot of the 
bistatic scattering versus 9 can be found. Each of these Figures shows the data for 
edge-on and face-on incidence on the left and right respectively, with the intensity on 
top and the phase below. The phase is that of a scattered field component at a large 
(constant) distance from the axis of the cylinder and is shown relative to that of a 
line source on the axis. The intensity plotted is actually the dimensionless quantity 
aie)/\ in dB. This particular normalization is convenient for computation and pre- 
sentation purposes, but since a\0) is a very complicated function of X (through, for 
example, the refractive index), it must be borne in mind that o{6)/X is not a wave- 
length-independent quantity. 

Spectral information is presented in Figs. 31 through 35 where o(0), o(7r), 

<y , <j and a are plotted (in dB ^.tm) versus X. Note that the explicit factor X 
T A ext 

has been removed, so tha' ‘-ere the 'normalization' is relative to a micron (qm). 

2 . 1 Angular Data 

The curves in Figs. 3 through 30 are self explanatory. They clearly show 
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Table 1: Computed Data for Hexagonal (Shell) Cylinder 
3 /im on side, 0. 1 /xm thick 


Xium) Pol. n Ret R/Z, 


0.76 

E 

1.307 

IP 

il. 70784 

3 

H 




4 

1.61 

E 

1.293 + iO. 000365 

IP 

0.00536 + i3. 81395 

5 

H 




6 

2.0 

E 

1.291 + iO. 00161 

IP 

0.02977 + 14.77438 



H 





2.25 

E 

1.254 + iO.OOl 

SW 

0,02740 + i6. 25466 

7 


H 




8 


E 

1.278 + iO. 000213 

IP 

0.00486 + 15.65462 



H 





2.6 

E 

1.206 + iO. 00080 

IP 

0.03876+ 19,10650 

9 


H 




10 

2.8 

E 

1.152 + 10.0123 

IP 

13.5260 + il. 17250 

11 


H 




12 

3.0 

E 

1.130 + iO.2273 

IP 

7.79740 + 13.41910 

13 


H 




14 


E 

1.045 + iO.429 

SW 

5.26973 - iO. 54082 

15 


H 




16 

3.1 

E 

1.280 + iO.3252 

IP 

4.20510 + 12.69040 

17 


H 




18 

3.3 

E 

1.530 + i0.0625 

IP 

0.55090 + 13.84930 

19 


H 




20 

3.5 

E 

1.422 + iO.0163 

IP 

0.24680 + 15.44060 

21 


H 




22 

8.0 

E 

1.312 + i0.045 

SW 

2.82939 + 117.23610 

23 


H 




24 

9.0 

E 

1.269 + i0.043 

SW 

4.09006 + 122.80586 

25 


H 




26 

11.0 

E 

1.093 + iO.242 

SW 

31.03981 + 17.98514 

27 

H 




28 

12.5 

E 

1.387 + iO.422 

SW 

12.08921 + 17.70071 

29 

H 




30 
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how the number of maxima and minima in 0 < 0 < y increases with decreasing X. 
and at wavelengths which are much longer than the side length of the hexagon, i.e. 

\ ^ 8 jum, the cross section has almost no angular structure. Changing the incidence 
from edge-on to face on has most effect in 'the directions closest to backscattering, 
and we also note the substantial differences between the results for E and H polari- 
zations. 

* 

2.2 Spectral Data 

The particular cross sections o(0), a{ic), o-_, cr^ and o are shown as 

T A ext 

functions of X in the infrared range in Pigs. 31 through 35. Each Figure has two 
parts, covering the ranges 0. 76 to 3.5 jam and 8 to 12.5 i^m. Separate curves have 
been included wherever the results for ec^e-on and face-on incidence are clearly 
distinguishable. lor the most part this is only with the backs cattering cross sec- 
tion or for X ^ 1 . 5 /um, and in other cases the differences are eonfiaed to the imme- 
diate vicinity of local maxima or minima in the data. 

In the shorter wavelength range most of the data were computed using the 
IP values for the refractive index. This range covers the main absorption band 
centered on X = 3 /jm and a secondary one at X = 2 jum as seen in the IP data 
plotted in Fig. 2. The wavelengths close to these show the main discrepancies 
between the IP and SW data, and at X = 2.25 and 3jum we theiefore ran compu- 
tations using both sets. The different refractive indices at X - 2.25 ,um do indeed 
produce substantial differences in the cross sections, and because of this sensiti- 
vity, we have not extended our detailed computations through the third absorption 
band centered at X ~ 1. 25 /um in the IP data. For X < 1.61 /um the scattering 
has been computed only at the single wavelei^h X = 0.76 jum of interest for a 
particular application. Since the IP data are given only for X 0. 95 )um, the neces- 
sary refractive index was obtained by extrapolation. 

The SW data for the refractive index were used in the longer wavelength range 
8.0 < X < 12.5 jam. 

The geometrical effects are particularly pronounced for 0.76 < X < 3.5 pm. 
This is not surprising since the dimensions of the cylinder are then comparable to 
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the wavelength, or are low multiples thereof, and this is the region where resonance 
effects and other interactions between the various contributors to the scattering are 
most important. As an example, while has a strong local maximum near the 
maximum in n. at >. = 3. 075 nm, the shape and overall width of the maximum in 
a is apparently affected by the side length of the hexagon being close to this wave- 
length. a(0) and a( 7 r) both show a corresponding drop in this absorption region. 

The behavior is quite different near the secondaiy maximum in at X = 2 /jm. 

For H-polarization but not for E, is large as expected, while o{0) and oiit) have 
local maxima for X just above 2 /uxo. with both polarizations. 

Further evidence for the way in which a geometrical effect can dominate a 
material absorption effect can be found by comparing the absorption cross sections 
at 2. 25 um computed from the IP and SW data. At this wavelength the SW value for 
n. is roughly five times the IP value, with n almost equal in both listings, but the 
SW value produced an absorption cross section which is about eleven times less 

than that given by the IP value for the refractive index. 

It is therefo'f'e obvious that any predictions of absorption and scattering 
based only on the properties of the material of which the scatterer is composed mav 

be considerably in error. 
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3. Conclusions and Suggestions for Future Work 

In this Report we have derived the theoretical basis for determining the scat- 
tering and absorption of electromagnetic radiation by thin-walled cylinders of infinite 
length and arbitrary cross sectional shape. Numerical procedures have been developed 
and have been used to obtain data for hexagonal cylinders which model sheath crystals 
of the type found in cirrus clouds. The procedures are economical for cylinders 
whose cross sections are not more than about 15 wavelengths in perimeter, with the 
cost decreasing rapidly with decreasing size. 

For wavelengths comparable to or less than the face length of the hexagon, 
the results (particularly for the back scattering) are quite sensitive to the polariza- 
tion and direction of the incident plane wave, i.e. on whether the field is incident 
edge-on or face-on as shown on p.6. At any given wavelength, the results can also 
be very sensitive to the refractive index employed Geometric effects can so influence 
the absorption that it is not at all safe to assume that the absorption versus wavelength 
curve will follow that of the imaginary part of the refractive index. 

The time available for the present study did not permit an adequate investiga- 
tion of the effect of wall thickness, nor allow us to do a detailed comparison of the 
results with those of the Mie-type series for hollov/ cylinders circular in cross sec- 
tion. Although our data for the bistatic scattering cross section versus angle are 
somewhat similar to those previously published (Liou, 1972a) for solid circular 
cylinders at near broadside incidence, data showing the precise role played by the 
geometry should have a high priority in any future continuation of the study. With only 
minor modifications our computer programs can also handle irregular additions to the 
hexagonal surface, thereby simulating rimed crystals, and allowing us to compute the 
effects of riming. Since the theory for non-broadside incidence has been derived, we 
would also like to develop the computer programs necessary to obtain numerica' data 
in this more general case, and once these types of data are in hand, it would be pos- 
sible and desirable to examine the forms of averaging that could sunplify the practical 
applicability of the da' a without losing its essential properties. 
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Fig, 31: Computed backscattering cross sections of hexagonal (shell) 
cylinders. In this and the following four figures, *and o 
show edge-on and face-on results respectively for IP data; 
similarly x and © are for SW data. 
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Fig. 32; Computed forward scattering cross sections of hexagonal (shell) 
cylinders. 
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Fig. 34; Computed absorption cross sections of hexagonal (shell) 
cylinders . 
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Appendijc: Scattering by a Cylindrical Resistive Shell 


A-i Formulation 

The shell crystal is simulated by a hexagonal cylindrical shell of uniform 
thickness and constitution, which is in turn treated as an infinitesimally thin, elec- 
trically resistive membrane. 

The concept of such a membrane arises nat'irally from a consideration of a 
thin sheet of highly conducting material whose permeability p is that of free space. 
If a is the conductivity and r is the thickness, we can define a surface resistivity 
R as 



ue.v T 
u e 


( 1 ) 


where x is the electric susceptibility, is the permittivity of free space and a 
time fector e“^^ has been assumed; and as f->0 we can imagine a to be increased 
in such a manner chat E is finite and non-zero in the limit. The result is an idealized 
(infinitesimally thin) electricafiy resistive sheet whose electromagneuc properues 
are specified by the single quantity E. In terms of the (complex) refractive index n 

of the layer material 


(n - l)kT 

where k and Z are the propagation constant and intrinsic impedance of free space, 
respectively. 

Mathematically at least, the membrane is simply an electric current sheet 
whose strength is related to the tangential electric field via the resistivity E ohms/ m , 
which may be a function of position. Since ju = there is no magnetic current 
present and 


S ,(e'^-e") = 0 

A”" — 


( 3 ) 
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where the affices + refer to the positive (upper or v)uter) and negative (lower or 
inner) sides of the sheet and fi is the outward normal to the positive side. If J is 
the total current 

(4) 

and from the definition of the surface resistivity 

n/^(a^E-> = -EJ . (5) 

The tangential components of the electric field are therefore continuous across the 
sheet, whereas the tangential components of the magnetic field have a jump discon- 
tinuity J directly related to the electric field via the resistivity R. With R speci- 
fied, the conditions (3) - (5) define a transition problem for the electromagnetic 
field and were first used by Levi-Civita (see Bateman, 1915) in studies of charges 
and currents close to a conducting sheet. 

A-2 Scattering by a Resistive Membrane 

It is convenient to start by considering a general resistive membrane in 
three space dimensions illuminated by an arbitrary electromagnetic field. We treat 
first an open sheet for which a representation of the scattered field was obtained 
by Knott and Senior ( 1974) and then examine the case of a closed resistive shell. 

If S is an open electrically resistive sheet and we surround it by a 
closed surface S^, the scattered field at any point r outside S^ can be written as 

E®(r) = -r ikZVj.jr"' 

( 6 ) 

H^(r) = V^VA/-ikYV,,jr 
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where Y = Z aud ir and are the electric and magnetic Hertz vectors defined 


= ^ \U^HgdS' 




v\EgdS' , 


8 == gtrlr') = 


with d = j £“ r'j is the free space Green function. On collapsing to the two sides 
of S, the expression for jr(r) reduces to 


T(r> * -Y \ -H")gdS' 


i. e. , x(r) = -^ \ \ J(rO g dS' 


where the integration is now a one-sided integration along S. Since 

^ - - -* — ^ 

nyi^(E -E") = 0, we have similarl3^ 

= 0 , (10 

and eqs. (9) and(10)are precisely those which would have been obtained fay starting 
with the concept of an electric current sheet. Hence 


E(r) = Elr)+V^V/,jr 


H(r) =H(r)-ikYV^ff 




where ^ , H are the incident field vectors. 
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If r is not on S, the derivative operations can be applied directly to the 
integrand of (9) giving 



In the particular ease v/hen the incident field and the surfece S are both independent 
of a Cartesian coordinate 2 , eq.{ 12) was the starting point for the analysis of resis- 
tive sheets by Knott and Senior <1974). Integral equations for the tangential conaponents 
of J were developed by taking the limit r— >S, and similar results can be obtained 
even in three dimensions . 

On the other hand, a somewhat different expression for the electric field 
is also possible. By applying standard vector identities to the second term in the 
integrand of( 12), we have (Senior, 1975) 



(14) 


where ie the surface divergence operator and the line integral is in the positive 
direction around the edge L of S. For a perfectly conducting surface the line integral 
vanishes by virtue of the edge condition. The electric lield expression ( 14) is then 
identical to that given by Poggio and Miller (1973) for a closed surface S, and since 


v'-J = ikYp = ikYn'-(E'*’-E’') , 
s “ s' 0 **■ — 


(15) 
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it is also identical to that produced by application of the appropriate Stratton- Chu 
formula (Stratton, 1941) to the closed surface consisting of the two sides of S. For 
a nonmetallic surface, however, the line integral is in general nonzero. The 
Stratton-Chu formulas are no longer valid and to avoid the occurrence of a line 
integral it is necessary to use the representation (12). 

If S is a closed resistive shell, the Stratton-Chu formulas can be used to 
give the following expressions for the total field at a point r outside S: 

E(r) = E^r) + ^ ^ |kz(n’^ H^)g + (n*/, e"^) ^ V’g + (S' • gj- dS' (16) 


H(r) = n\r) + \ \ ■[-iky(n'A e'*' ) g+(n'^ H‘^)AV'g+(n’ • )V'g| dS» ( 17) 


where the unit vector normal is directed into the space containing r. In contrast to 
(12) and (13), the surface integrals now involve the fields on one side of S alone, 
rather than their jump discontinuities across S. We recall, however, that (16) and 
(17) are obtained by application of the vector Green's theorem to the volume exterior 
to S containing the observation point. If, instead, we apply the theorem to the vol- 
ume interior to S with the observation point still outside S, g will be regular through- 
out the entire volume of integration, implying 



By choosing r the same as in eqs. ( 16) and (17), subtraction of ( 18) from (16) gives 




where |e] " E^ - E~ and = h"^ " H . Similarly 


52 


I 


1 


I 


rr 

H(r) = hV) + /i 


)g + (n*AiH])/v'^'g+{^* [H])V'g dS' 


and when the boundary conditions C3), (4) and the relation ( 15) are used, we have 

E(r) = l\r) + ^ ^ - -^(V^ • j)Vgj dS' (20) 

d J k 

H(r) = H^(r) + ^ ^ J(r’) V'g dS' . (21) 

The above results could also have been obtained using a Hertz vector repre- 
sentation of the scattered field, and are identical to those for an open resistive sheet 
when the line integral contribution to the electric field is ignored. More to the point, 
when the steps leading from (12) to (14) are reversed and applied to (20), eq. (12) 
for an open sheet is recovered, thereby validating this integral representation for 
both open and closed resistive sheets. 

For a closed shell either (12) or (20) can be used to generate an integral 
equation for J, but the two equations are significantly different. On selecting the 
tangential components of (20) and then allowing the observation point to lie on either 
side of S , we have 




J(r')g - -2 (V ■ J)6»V'g|dS' 


and since 

^AV'g = (n - a')A^'g+ 

the more singular term in the integrand contains only tangential derivations of the 
kernel at the self point r' ~ r . It follows that the integral in (22) is continuous as 
r approaches S, allowing us to apply the limit directly to the integrand, and impo- 
sition of the boundary condition (5) then gives 
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fiAE^r) = Ri\\J{r) - ikZ 



(23) 


for r on S. This is a valid integral equation whose only disadvantage is the occur- 
rence of surface deru'atives of J, 

This difficulty can be overconae if we use the integral representation ( 12) . 
On paralleling the steps leading from (20 to (22) , we have 




k 


(24) 


and because of the higher order, non-integrable singularity of the integrand at the 
self point, it is no longer possible to apply the limit directly to the integrand. Since 
the contribution of the 'self cell' tends to infinity as the cell size tends to zero, even 
the limit shown in (24) does not exist, and though (24) differs from (22) only through 
an integration by parts applied to the second term in the integrand, it does rot con- 
stitute an analytically valid integral equation. Nevertheless, it can be used as the 
basis of a numerical solution provided the segmentation of the integral which is 
inherent in such a method is performed prior to the limiting operation, with the 
segment size remaining non-zero. The resulting equacion is in some respects 
preferable to (22) and has been found more convenient for our present purposes. 


A-3 Scattering by a Cylindrical Resistive Shell 

We now turn to the problem of a closed cylindrical shell illuminated by a 
plane wave at oblique incidence, and consider the form which the above integral 
equations take in this particular case. 

Since the shell is independent of the z coordinate, the entire dependence on 
z is that produced by the incident field. If, therefore, 

E^r) = E\p)exp(ik z), H^(r) = H^(p)exp(ik z) (25) 

^ “ ““ z Z 
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where r — £ + zz, the total fields i.dmit the same decomposition and, in particular, 

J(r) = J(s)exp(ik^z) (26) 

where s is the circumferential distance around the shell in the plane z = 0. Thus 


where (now) 
and 


= -^expUk^z) J(s' 


)H^^^{K,d)ds' 


d 

K « . 

Z 


(27) 


The remaining integration is with respect to the circumferential distance a around 
the (closed) perimeter of the shell in the plane z = 0. Also 

7; . J = exp(ik^.-) J^(8’)+ik^y s')j 

and since 


^ J K£') exp(ik^z')7>g dS' = I exp(ik^z) ^ f(e’)(V» - ik^5) nj^^^Kd) ds‘ 

where ^ is the two-dimensional (transverse) del operator in the plane z = 0, 
eq. (23) becomes 

= Rca£(s) + j' |£(b’)h[j^W - -| 


r 9J 

+ ik J (s') 
ds' z z J 


•(V| - lk^S)Hjj^\xd) (. ds' 


The S component is simply 


Jc ft as' J ^ 


(28) 


(29) 
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which is a. coupled integral equation involving the longitudinal and circumferential 
components of the current. Likewise, from the z component of (28) we obtain the 
second coupled equation 


Y1e[{b) = YRJjs) + 7 \ ] JJs')(s . s')Hj.^^(fcxl) 


-4 777 + JJs') (s-W/ Vdl^ds’ . ( 30 ) 

j^ 2 {^ 8 S' Z Z j 1 j 


and though it is not in general possible to decouple (29) and (30), we can eliminate 

J from the integral portion of (30). Differentiating (29) with respect to s , we 
z 

have 


J (8»)(^^)H;^\Kd)ds' = 
z 1 3 OS 




Jc 



fc r aj ... 

- i“| \ (s *d)H^ (Kd)ds’ 

Jc 




and when this is substituted into (30), the integral equation becomes 


. K 




RJ (s)-i4 — RJ (s) 
8 2 ds z 

K 



dJ 


fj (a')(s-s')Hj,^Vd)--^ ^(S 

S 0 K ds' 


A ( 1 ) 

d)HY^(Kd) 


ds' 


( 31 ) 


which now involves only as a pseudo excitation terra. 

Formally at least, the corresponding equations resulting from (24) can be 

obtained from (20) and (30) by integrating by parts the terms involving 9J /9s'. 

z 

Equation (29) is then replaced by 
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YE^s) = YRJ (s) +~ ] fj (s')Hl^^(/cd) - i— J (s')(i* • &)Hl^Vd)] dS' (32) 

2 Z 4bl£ Z U /C S J. J 

c 

and (30) by 

Ye\s) = YRJ (s)+x \ J (s')(s-s')Hf,^^(Kd)ds' 

s s 4 ^ s 0 

C 



The latter is meaningful only if the second integral on the right hand side is handled 
in the manner described earlier. Wdth this proviso, however, (32) and (33) are 
adequate for a numerical determination of the currents and we can, of course, elimi- 
nate from the integral portion of (33) if we so desire. 

The case considered now is that in which the plane wave is incident in a 
plane perpendicular to the z axis of the cylindrical shell. Th^n k =0 implying 

X = k, and (32) and (33) reduce to 


YE (s) = 


YRJ (s) +- 
z 4 


J (s')H5,^\kd)ds' 
z u 


(34) 


YE (s) 
s 


YEJ (s) +7 
s 4 


J (s')(s.a')Hi^\kd)ds' 
s 0 


1 fim 

+ 4£,-^C 


j' [(s-5)H^^^\kd)] ds' , 


( 35 ) 
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which are two uncoupled integral equations for the longitudinal and transverse current 

components J (s) and J (s) respectively, 
z s 


If the incident field is E - polarized, i.e. 

gi _ A^-ik(xcos0Q + ysin0^) 

where 0^ is the angle of incidence with respect to the negative x axis (see Fig, A-1), 

only the component J (s) is excited and we shall refer to (34) as the E - polarized 

z 

equation. O^ce J (s) has been found, the scattered electric field in the far zone is 
z 


e“ = 2 (e 9 ) 

'/rko e' ' 


/rkp 'E'“-’''0' 

where the complex scattering amplitude is given by 

E 




J (s')e 
z 


-ik^ • £' 


ds' 


(36) 


with 


p = xcos 0 + ;^sin0 


and in terms of P, the two dimensional scattering cross section is 

a(e,So) (37) 

If, on the other hand, the incident field is H- polarized, i.e. 

„i _ A -ik(xcos0 +ysin0J 

Jti “ 2G U (J 

the only component excited is and this can be obtained from the H - polarized 
equation (35) . In the far zone 


rrS A [2 i(kp - 5 t/4) ^ , 

If = zJ— e ' P„0,0J 

^ TTkp H ' 0 


where 




k 

4 


J' p • n'J^(s')e ^ ds' , 


(38) 


from which the H- polarized cross section can be found using (37). 

Two quantities of particular interest are the total (integrated) scattering and 
the absorption cross sections 0 ^ and respectively. The former is given by 
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' 2; 


i2jr 

(7(0,0g)d0 

'0 


(39) 


and from the forward scattering theorem the absorption cross section is 

= - a^(6g) - I Re. P(0g + 7^,6^). (40) 

We remark that for a finite shell of length £ (»A,) each three dimensional cross 

(3) 

section o for incidence in a plane perpendicular to the length and computed on the 
assumption that the surface field is the same as for the infinite shell is related 
to the corresponding two dimensional cross section by 



where X is the free space wavelength. 


(41) 


A-4 Computer Programs 

The integral equations (34) and (35) are special cases of those solved by 
computer program RAMVS (Liepa et al, , 1974), and because of the availability 
of this general program we chose to concentrate on these equations rather than 
the ones involving the derivatives of the currents. While seeking to refine RAMVS 
and to make it more efficient for a polygonal shell, we became aware of certain 
errors and/or deficiencies in the program which are most apparent when the 
resistivity is small. These had not shown up in the testing and verification done 
earlier, and it proved quite time consuming to locate the errors and rectify themi 
We believe this has now been done. In addition, the program has been extended 
to compute the normalized two dimensional cross sections and to plot the normalized 
bistatic scattering cross section, where the normalization is with respect to the 
free space wavelength to make the quantities dimensionless. The two programs 
that resulted are designated RICE and RICH and are based on eqs. (34) and (35) 
for E- and H- polarization respectively. 
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Both equations are solved by breaking up the integrals into N equal segments or 
cells within each of which the current is assumed constant, but because the kernels 
are infinite when the integration and observation points coincide, the self cells 
must be treated anal}rtically. This is a rather trivial matter in the case of eq. (34) 
whose kernel has an integrable singularity. Prom the small argument expansion of the 
the Hankel function we have 


.( 1 ) 


2i 


kd 


H^^-^^{kd)^l + ^(in^+7) 


0 


7T 




2i 

r 

7r(kd)‘' 


where 7 == 0.577215. . . is Euler’s constant, and hence, for the self cell A of width 
26, 


J^(s')HQ^\kd)ds’ Ci 2&r^(s) [ 1 + ~ [fn 0. 028798. . ,J j 


Equation (34) now becomes 

YE \s) = [ye + Y ^ Y + 0.028798. . . ] j J 


,(s) 


(42> 


+ JL 
2X 


J (s')Hl^\kd)ds’ 
z 0 


>C-A 


and this is the equation used in program RICE. 

The same reasoning applied to the first integral in (35) gives 

J (s')(s.^’)H[,^V<d) ds’=:^26J (s) fl+ — (in ^ + 0.028798... 
s 0 S L 7T ' X 

A 

but the second integral is more difficult to treat. When the self cell is excluded 

from the range of integration, the limiting operation can be applied to the integrand 

directly, and since J (s’) is assumed constant over each cell, 

s 
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iim 

£-*C 


Jg(s')|^,[{S-a)H^l^Vcd)]cls- = 


N 


X-A 


P’ =P, + 6 


^ Jg(s ) (s-^)E^^\kd) 1 ^ 

j = i -^^'=0.- 

j/i 


p'=P-6 


where the subscripts j and i denote the integration and observation points respec- 
tively. To evaluate the self cell contribution we first use the fact that 

~[(s-W^^W)] = k[(n'.a)(n.&) - {s'.S)(s.d)}H^^^’(kd) 

- k(n'*S){n-a)Hjj^\kd) 

with the prime attached to the Hankel function denoting differentiation to express the 
contribution as 


where 




I = k 


^ ' (n'*&)(n-d 


d) - (s'-d)(s.d)j H^^^^'(kd)ds'. 


For a locally plane element with the observation point a small distance y above 
its midpoint, 

I 

I = 2k 


'pQ y + t 


and on inserting the small argument expansion of the Hankel function, we have 

6 


—L_ 

TTk 2 ^ 2 

y +6 


2 .2 


Z. 




dt • 


\ 
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Hence 



0 


the first term of which becomes infinite as the cell size tends to zero. Equation (35) 
now has the form 

YB>) = U^-|-|j+||{lff(Aif+ 0.028798...-0.5)}]j,(s) 

*• 7T 



whose solution is computed by program RICH. 

Of the two programs it is evident that RICH is the more complicated, and 
to obtain a feeling for the rapidity of convergence as a function of the cell size, 
the program was run for edge-on incidence on a hexagonal shell 3p on the side at a 
wavelength 2 . 6^. The normalised resistivity employed was 

YR = 0.03876 + i9. 10650 

corresponding to a shell 0. Iq thick and the refractive index 

n = 1.206 + iO. 00080 

for ice at 2. 6ju quoted by Irvine and Pollack (1968). The results of decreasing the 
total number N of cells from 72 down to 24 are listed in Table A-1, and we observe 
that the total scattering and absorption cross sections are relatively insensitive to 
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cells per 
wavelength 

time used 
(sec. ) 

(dB) 

a<5r-H0 0)/X 
, 0' 0 
(dB) 


■Ja(9o>A 

10.4 

10.911 

“44. 83 

-6.45 

0. 02229 

0. 00158 

8.7 

6.645 

-44. 70 

6.45 

0. 02230 

0. 00157 

6.9 

3.768 

-44. 48 

-6.45 

0. 02233 

0. 00155 

5.2 

1.948 

”44. 04 

-6.46 

0. 02241 

0. 00153 

4.3 

1.332 

-43. 63 

-6.46 

i 

0. 02250 

0. 00151 

3.5 

0.895 

-42.91 

-6.46 

0. 02272 

0. 00149 

















n r> n o 


N, as is the forward scattering cross section, but the backscattering changes 
noticeably on. decreasing the sampling rate below (about) 9 cells per wavelength. 

In running the program we have therefore used 9 or more cells per v/avelength 
whenever possible consistent with the maximum of N = 100 aliov/ed by the matrix 
Inversion routine used. . • * 

Listings of the two programs follow. 




• «!,. aJ* aAa aAa ala ^ 

«a,a a,aap..a,» a^a a,m a,a af» ay» aT* a^r* aja a»a aja ag. 

INPUT FfIRMAT FOR RRH^HAM RICF SFPT, 1975 

tsfs »J* Uaala Uaala-afci aJa aXia>a tJa aJa al. .laat. aXa aU >laaJa a»a a>. JLa afa^la aW ala ala ^ aU^aia .la ala a»a aia aJa alia .jC. - 

'•T* **• I’* »r» •V' ^ » t < i *r» *t* 'f ■ «*t^ ,^r» .j" •»!» * 

CARn 1 FnRMAT {1RA4) TITLE CARD: USE UP TH 72 CnLUNMS 


C 
C 

C CARD 2 
C 
C 
C 

c 
c 
c 
c 
r. 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


FHRMAT ( I?, I3t 5Fin.5J MpRFtKnnF, ZFAC , UAUe, FIRST, lASTt IN< 

MpRF=n THIS WILL PF THF LAST RUM FUR THIS DATA SFT 

MORFrrl THFRF AR F MRR f pATA TP RE RFfip AFTFR THIS SFT 

KPPF=n GP-MPIITFS RTSTATTC .scatter I.’"P pattfrm 

KnnF=i cnMPHTFS rackscattfr iwc pattfrm 

ZFAC A COMPLEX FACTOR MULT TPLY IMP ALL ELFMENT 

WA\/F WAVFLENPTH 

FIRST INITIAL SCATTFRIMP AND INCinFMCE ANGLE 

LAST FINAL AMPLF 

INK AMPUL AR TNCRFMFNT 

CARO 3 FORMAT t I 2 , 5 X , 7F IP .5 } N , Z . X A , Y A , XR , YR , AMP 


CARO 


N 
Z 

XA» YA, XRtYR 
AMP 


mijmrfr pf samplimp PPIMTS pn this segment 

MPRMALIZFH IMDFDAMCE Or THF SEGREMT 
SFGMFMT FMOPniMT.S 
angle SURTEMOEO PY the SEGMFMT 
ZERO IN COL ? SHUTS OFF 
REAOIMG PF segments 


^ EPRMAT ( I2» I3aFin.5 ) MORF , KPPF , ZFAC tF 1 RST rL AST , INK 

THIS CARP TS US EO only IF, 

PM CARO 2, MPRF=I 


in SPECTFICATIOMS: 

5:INPIIT DATA: A:nilTPUT{ PR INTER ) ; 7 : nUTPUTf GR APH ICS ) 

f' ?r» ^aia ma2s!it2l!£s!£ tissis 2ft ^ jJj *2j 2 *f *^ *•* **^ “** ^“*T 2 !t 5 ?f *^‘* *^ “** *«‘» ■^•^* '•'•♦ta- aJa*kU ataaJa Oa ala ala ala ala. U ala ala ala ala a Waia .JU ala ala aJa al*. ala ala ala ala «la a'aala aU ala ala ala al. 

\_g *T»*V •!— •'■"•T’ 'r?rr* a^^aga.agaaf.a^agar^anaay a^a.agaa|a.a^aaya agv a^a 2^2^ a{C^a aga a«ia^ a^ aft ajt ^ 


w&'f-ieBis 
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o o P» o r) n o n o o o n o o n n o o n o o r> n n n o o n n o o o p 


GnMPLFX=!=a A( I 0 a,i 01 ),PHI( 10 O)TPIMKC 100 )»ZS{lO) 

CnMPLFX-R DFLtSUI^TZFAC 
RFAL-4 LAST* IMK 

ReAL?!:A X<-2nO)*Y( 20 n) »XM{200) »YM( 20n) ,S(200)TnSPl 10) tAStJHOAl) 
TMTFGFR=!^4 ID( IR) ♦ LUMP ( 200) ♦ TPnL(? ), LV( 100 ) ,MM( 100) 
data MH/ion/ 

CnMMnN YPIFS/ HI»TPI, PIT, PIPIt YZ,REn,nlG 

DATA. IPDL /'FFFFVr 
DATA IPP/1/ 

C....-RFAD TMPiiT HATA AMD GFMERATF RPDY PROFILE . 

5 REAo ( 5, 100) rn 

rFAD ( 5 , 200 ) MnRF,KnnE,ZFAC,WAVE>FIRST»LASTT INK 

wRiTE(7,ioo) in- 

WRITEITtIOI) i-'AVF 

WRITE(7,19P) FTRST,LAST, IMK,IPP 
iOl FORMAT! ' LAMRnA=» ,F5. 2,' MICROf^S* ) 

199 FORMATC lXt3Fin,5, 17) 

I WHICH=1 

rF(REAL(ZFAC),FO.O.AWn,ArMAO(ZFAC).EO,n, ) . 

K Z FAC= ( 1 . F -2 n , 1 . .t -20 ) 

writf ( 6, 150) in 

WRITE (6,300) 

CALL OEnM(Ll!(-^F*X,Y,XM-,YNTSTnSOTZSTM) . ■ 

MT=M/2 • 

LL=LUMP(M) 

20 IF, (KnnF.NE.O) 00 TO 25 

MTMC=1 

MRIT=1+IFIX( (LAST-FTRRT)/IMK) 

on TO 30 • 

25 . WB;! T=0- ■ ' 

WTMC= 1 +I FI X( ( LAST-FIRST) /INK) 

30 write { 6,1501 m . 

• WRITE (6,400) IP0L( TPD),LL,MT,NTMC,MRIT,OAVE,ZFAC 

WRITE (6,4^251 first 
. 38 on 3 5 1 = 1 1 L L 

ZSM)=ZS{ I )5=Z-AC 
35 nSD( I )=nSQ( I ) /wavf 

XK = TPIZWAVF^ 

c .... .comstruct matrix ELFMFNTS 
IF(IWHlCH.Fn. 2 ) on TO 37 

CALL MTX (M,MW,XK..X.Y,XW-yri,n.SQ, LIIMP,ZS, A) 

C.,...COMPUTF iMCinFMT fIFLO AMD INU/ERT MATRIX 
37 TFTA=RFni!sFIRST. : 

CT=CnS( TETA) 

ST=SIM(TFTA) 
no 60 1=2, M, 2 

HOLn=-XK*{CT^'X( I )+ST=i=YC I ) ) 
nFL=CMPLX(COS( WOLny, STM (NOLO) ) 

60 PINK { 1/2 )=r)FL 

CALL FLIP( A,MT,MH,LV.MM,PIMK,PH1 ,IWHICH) ^GE Ig 

“OE Poor QUipjjY 


RSUM061) 


I 


I 

I 


1 


1 





C PRIMT ntlT CURRENTS AND ELEMENT PROPERTIES FOR FIRST ANGLE 

WRITE (ft*50n) 

IT=0 

no 65 I=2»M,2 
IT=TT+1 

AMP=CARSI PHK IT) ) 

PHASEs 5DIG=:^ATAN2{ AIMAG(PHl(lT))fRFAL(PHT{IT))) 

65 write**! 6»250) ITtI5FG,X{I)tY{I)tS(T)»nSO(ISEG),ZS{TSEG) » AMP » PHASE 
C....«nnpE OUT THF APPROPRIATE FIELD FACTORS 
thf=first-ink 
K=0 

IF (KOnE.EO.I ) GO TO 70 
WRITE !6,aon) FIRST 
GO TO 75 

70 WRITE (6» 600) 

75 THE = THF+INK 

IF ( THF. GT. LAST) GO TO 105 

IF ( THE, EO. FIRST) GO JO R5 , ^ ^ 

C...-»IN THF FOLLOWING LOOP» PINK IS NOT NECESSARILY THE INCIDENT FIELD 

TETA=REn*THF 
CT=CnS{ TFTA) 

ST=SIM( TFTA) 
no RO J=2»M,2 

HOLn=-XK*(CT*X( J)+ST*Y{ J) ) 

DFL=CMPLX(CnS( HOLD) ,SIN(HOLO) ) 
pn PINK! J/2 )=DFL 

IF iKODF.EO.O) GO |0 R5 

CALL FLI P( AtMT»MH. LV, MM, PIMK» PHI » 2 ) 

R5 SIIM=CMPLX!0.0,n.O) 

C ADD UP THE CURRENTS 

DO 95 J=2,M,2 
JT=J/2 

95 SUM-SIIM+PHI ( JT)*PINK! JT)*DS0(LUMP! J) ) 

SUM=-Sl)M 
SUMR=RFAL! SUM) 

SIIMI=ATMAG! SUM) 

SI IM SO = SUMR''l'- SI IM R + SUM I SUM I 
SCAT = in.=i=ALnGlO! SUMSO ) + 1.9 612 
K=K+1 

ASUM!K )=SIJMS0 
RSUH! K)=REAL( SUM) 

PHA5E=ni G=!*ATAN2 { SMMI < SUMR ) 

WRITE! 6,901) THE , SCAT, PHASF 
WRITE(7,900) THE, SCAT, PHASE 
GD TO 75 

105 DIFF=LAST-FIRST ^ ^ 

IFiDIFF ,NF. IRO.O .AND. DIFF .ME. 360.0) GO TO 205 

inDD=H0n!K,2) ^ 

IF!DIFF .EO. 360.0 .AMD. lODD .NE. 1) GO TO 205 
IF! DIFF .ED, IRO.O) FAC-2. 0 
IF(DIFF .EO. 360.0) FAC=l.O 
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Cj!:*# SIMPSON IMTEGRATION OF CROSSFCTIONS 
KLAST=K 

iFd'nnn ,eo. o) klast=klast-3 

SIGT=ASIIM{ 1)+A5UM( KLAST) 

siG=o.n 

DO 203 1=2, KLAST, 2 
203 SIG=STO+ASnM( I ) 

SIGT=SIGT+4.n*SI0 

SIG=0.0 

KLAST=K-2 

no 210 1 = 3, KLAST, 2 
210 SIG=SIG+ASIIM( I) 

SIGT=SIGT+2.n*SIG 

SIGT=SIGT*IMK/3.n 

IFdODD ,E0. 0) SIGT=SlGT+3./R.=!=INK=:" 

S{ASllM(K-3)+ASnM{K) + 3.=:=C ASUM(K-2}+ASIJM(K-n ) ) 

SIGT=SI GT=!^FAC4REn/4. 
gMH SIMPSON INTFGRATIOM 
KF=K 

IFIFAC .FO. 1) KF=(K+l)/2 
SIGF=-RSIIM( KF ) 

SIGA=-SIGT4-SIGF 
IFiKOnF) 220,220,240 
220 SIGTnR=10.0*ALnG10( SIGT ) 

IF(SIGA ,LF. O.ni GO TO 225 
SIGAOB=10,0*ALnG10( SIGA) 

GO TO 230 

225 WRTTEt f.,825) S TGT, S IGA ,S IGF 
WRITE(7,R25) SIGT,SIGA,SIGF 
GO TO 205 

230 SIGTDR = in,0=^ALOGlO{ SIGT ) 

WRITE! ft,R50) STGT,SIGTOR,STGA,SIGADR 
WRITE(7,«51 ) SIGT, SIGTnR,SIGA,SIGAnR 
GO TO 205 

240 SIGTnR = in.n*ALriGlO( SIGT) 

WRITE{6,R75) SIGT.SIGTOR 

825 FORMAT! ///, 5X,40H4FnHL'!‘ WFGATIVE ABSORPTION CROSSECT I dW, / /, 
S15X, 5HSTGT=, FR. 5.5X .5HSIGA=,FR .5 .5X, 5HSIGF=, FR.5 ) 

850 FORMAT! ///, 15X, 5HS I GT= , F8. 5 , FB. 2, 3H OB , 5X , 5HS IGA= , F8. 5 , FR . 2 , 
&3H OR) 

851 F0RMAT(///,i5X, 'SIGMA(T)/LAMRnA=»,F8.5,FR,2,' ORi^SX, 
fi»SIGMA( A) /LAMRnA=» ,Fa,5,FB.2, • OR* ) 

875 FORMAT! ///,15X,7HSIGAVe=,FR.5,F8.2,3H OR) 

205 IF IMORF.EO.O) GO TO 5 
no 103 1=1, LL 
ZSd)=7.S!I)/ZFAC 

103 nsod )=nso! n^wAVE 
nEL=ZFAC 
PWAVE=WAVE 


PAGP rcj 
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READ (5*200) MORE , KOOE * ZFAC * WAVE * F IRST , LAST , IMK 

WRTTE(7*10n) m 

WRITE(7*101) WAVE 

urITE(7* 199) FIRST* LAST, IMK* I PP 

IWHICH = 1- 

IF{REAL(ZFAC),FO.n. AWn,AIMAG(ZFAC).Ea-0. ) 
flZFAC=( l.F-?0,l.F-?0) 

IF(REAL(ZFAC),FO.RFAL(nEL),AMn.AIMAf;( ZFAC ) . EO. AIP.AG{ PEL ) 
R.AMn.WAWF,Fa.P'''AVE) IWHICH = ? 
on Tn 20 

100 FORMAT (1«A4) 

ISO FORMAT (1H1,1RA4) 

200 FORMAT ( 1 2 , I 3 , FFIO .5 ) 

250 FORMAT {2I5*7F10.5*Fin,31 

300 FORMAT ( lOHOSEO IM, 1 IX *2-^OFM0POlMTS OF THF SEGMENT* 19X, 

S1RH5EGMENT P AR AMFTFR S / 1 IH 'M!M CF LL S * A X , 2 HX A * RX , 2HY A* RX , 2HX B , 8K , 
r.2HYR* AX*24HAM0LP RADIUS L FNGTH,4X * 14HRE-Z IM-Z/) 

400 FORMAT ( //?1 X , lAHKFY PARAMETERS// 

£16X*21hIMC10FMT POL aR I 7 AT T0M*?2X., 1 Ai/ 
fil<=.X,23HMUMaFR OF SFOMpMTS IISFO*!?!/ 

£1 AX*33HTnTAL OF CFLLS ON THF RnoYtlll/ 

£l^>X,35H/JliMRFR OF JMCIOFMT FIELD DIR FCT IONS * 19/ 

£1AX,29HMMMRFR OF FISTATIC 01 R FCT'T OMS * 1 1 5 / 

£1 6X, IOHUAVELFMGTH, F3A.. 5 / 

£1^X*1 EHIMPEOAMCF FATnR*Flf.,5* » ,».Fin.5) 

425 FORMAT (////, 3nX, 1 R^SUR FACr FIFLO nATA/,?7X, 

£29HF0R TMCIUFMT FIFtn 0 T R FCT I 0M= * F7 .? ) 

500 FORMAT (UNO .1 S-G-4X.AHX( I )* AX,-iHY( I )*5X,-:^HS( I ) *5X*ftH05Q( r 5 , 
£4X.FHRS(I) ,^-xi5HXS(I) ,4-X , ftHMHO ( .1 ) ,AX . fSHARO ( J ) / ) 

600 format ( lHl,27X*2P-3AC'<SCATTFRlf''0 CROSS S FCT I PN//23X * 

£36HTHFTA SIGM A/LAMBOA .OR PHASF*PEG} 

flOO format { 1H1*23X, 33HBISTATir SCATTFRTNG CROSS SFCTIDM/23X, 

£29HF0R IMCrOFMT FIFLO OT R FCT inM= , F7. 2 / /? 3 X * 

£36HTHFTA S IGM A / L AWFO A , OR PHASF*OFG) 

900 FORMAT ( 15X,F13.2*F14.2*F16.1) 

901 FORMAT ( 15X*F13.2*F1A.,2*F15.1 } 

FMO 

SUBROUTINE MT X ( M , MH , X K . X, Y * XM * YM, OSO, LUMP * Z5* A ) 

C RICE VERSION, --POLA ' I 7.ATI0M 

RFAL*^- X( 1) *Y( I ) , '{ 1 ) , Y^M 1 ) ,DSO( 1 ) 

INTEGER*^*- iJiMP(l) 

COMPLF X=i-fi A ( MH, 1 ) , Z S ( 1 } * W2 , ni i.M 

common /pies/ PIfTPI, PIT* PIPT, YZ*REn*niG 

IH=0 

no 10 I=2*M*2 

IH=IH+1 

XI=X(I) 

Y!=Y(I) 

JH=0 

no 10 .1=2, M, 2 
JH=JH+1 
ISFG=LIIMP{ J) 
nS=nSO(ISEG) 

IFnH* EQo JH)GO TO 50 
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Rx=xi-x(.n 

rY=YI-Y(J) 

R.?SORT{ RX'^RX+RY’i'-RY) 

R K sR K 

CALL HAiVKZK RK ,n,HZ.nilM) 

IH,JH )=PIT=^HZ-nS 

Rn in 10 

50 A( IHt,IH)=ZS( ISFG) + 

£DS*CMPLX( PITtAinOt OS ) + .02B7C3R37 ) 

10 CPMTINIIE 
RFT1IRW 
EMO 

SIIPROUTIME RFHM { LUMP » X . Y ,XM ,YM.S f nSQ«ZS *M ) 

C THIS VERSTHN RFAHS ANT) RFMETATES SEGME^'TS IN CDUNJER-CLOC KW I S E DIRECiIDN, 
C THE SURFACE MUST RR CLUSPn. 

C pniMTS ARE RFMFRFTPn AT THE START AMD MIOPDIMTS OF EACH CELL; 

C THE START POINT OF THE FIRST CELL EVENTUALLY CPINCIDES 
C WITH THE FNO POINT OF THF LAST CELL. 

CnMPLEX*B ZS(l).-7 : 

X ( 1 ) . YC 1 UXM( 1 ) . YMf 1 ) ,S( 1 ] .PSO( 1 ) 

TMTEGER^A- LOHP(l) 

COMMON /PIES/ PI^TPI-, PIT. PIPI. YZ.RPn.nlG 
1=0 
L=0 

C.....R^AD INPUT parameters AMD PREPARE TO GENERATE SAMPLING POlhiTS 
’’io'rFAO (S.200) N.ZtXA.YA.XR.YR.ANG 
IF(N.LT.I) go to 120 
LIM=»2*N 
TX=XA-XR 
TY=YA-YR 

n=SORT( TX=!=TX+TY*TY) 


L = L + 1 
ZS{L)=Z 

IF (AMG.FO.n.O) on TO 20 

T = n , 5 * RFn=:sANG 

TRX = TX + TY*CnT Af!f T ) 

TRY = TY - TX 5 :: C 0 TAM ( T ) 
RAp = 0 .‘?* n / SIN(Tl 
ARC=2.0*RAO-T 
ALE=T/N 

nso( U= 2 .o*rap=:'Alf 
GO TO 30 
20 RAn=9S9.S99 

ARC=D 

nso(L)=n/N 

C..., .start GENERATING 
30 nss=o. 

ns=nso( L )/ 2 . 

no 50 JJ= 1 .LIM 
1 = 1+1 
J=JJ “1 
sm=pss 
nss=DSS+ns 


OBIGWAK PAGE IS 
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LlJMPd ) = L 

IF (/SNR, FO. 0.0) on TO 40 

ARR=J’.’=ALF 

5IMO=5IN{ ARR) 

CnSQ=CnS ( ARC? ) 

X( I ) = XA-n.5* [TRX#( 1 . 0 -r.nS 0 )-TRY*SIND) 

Y( n=YA+n.5-4(TRX*SlN0+TRY*( l.O-CnSO) ) 

XN(I ) =+0.5*( TRX^^.n50+TRY*RT^'0) /RAO 
YM( I )= n, 5* ( TRX*S [MO-TRY-CnSO ) /RAO 
Rn Tn 50 

40 X( I )=XA-0,5*.I*TX/M 
Y( I )=YA-n.5*J-TY/N 
XM(I )=-TY/n 
YN(I)= TX/n 
50 COMTINUF 

ion WRITF( A, 300) L,NtXA,YA,XR ,YR , ANR.RAn, ARC » 7 . 

an TO 10 

120 M=I 
LL=L 

200 FriRMAK I2,3X.7Fin.5) 

300 FnRf-'.AT( 13, I 6, 4X , 4F10. 5 , F7 . 2 ,4F1 0. 5 ) 

RFTIIRN 

ENO 

SIIRROIITPJE HAMKZl ( R.N.HZERH.HnME ) 

HANKEL FliNCTinHS ASF OF FIRST KIMO — J + IY 

...a. N = n RFTltRMS H7FRH 
. .... 't = l RFT!!Rf'S = 

M=2 RPTORNS H7FR0 AWn WOMF 

SliRROnir'E RFniijRFS R>o 

SI!RR^HTI^iE AOAM f‘‘ll5T RF SllPPLIFP RY USER 

D I MENS I dm a ( 7 ) , R ( 7 ) , C ( 7 ) ^ n ( 7 ) t F < 7 ) , F ( 7 ) « F ( 7 ) , H C 7 ) 

COMPLEX HZFRO.Hnric 

DATA A, R, C, D. F. F, r,,H/l.n,-2.74og9q7, X ,;?F-5A20a, -0. 31A3RFA, 

SO ,0444470 , -0 .003O4A4 ,0 ,0002 1 ,0 .3R74(^F,Q 1 , 0. F055O3RFI , -0. 74350384., 

SO. 253001 1 7, -Fi, 04 7 ^-1214, 0. 0ni-27Oln,-0, 00024P.4^, 0. 5 ,-0. 5F.2499R5 , 

FxO, 21 00 35 73, -0.030542 RF , O , 0044331° , -O .00 031751.-0.0000110°, 
f.-O. 53 5519 8, 0. 2212091 , 2. 1 58 270°, -1 . ? 15482 7, 0. 3123°5 1 ,“0. 0^00075, 
8*0.002 78 73, O. 797.5 845 <^,-n .00000077,-0.0055274, -0,00009512. 

SO. 00137237, -o, 0007 2P 0 5, 0. 00014475,-0, 785 79 816,-0,04 155397, 
8.-0.000039 54 ,0.00 2 57 5 73,-0,00 054125 ,-0.000293 33,0.00013 55 8, 

SO, 797884 5 6, 0. 0000015 R, O , 0 1 5 5° 557, 0, 000 1 7 1 05 , -0, 0024 95 1 1 , 

SO. 001 13 553 ,-0.000200 33,-2 .35 61°‘949 ,0.124995 12, 0.0000 565, 
S-0.00537879, 0.00074348,0,00079824,-0.00029156/ 

IF (R.LF.O.O) 00 TO 50 
IF (N.LT.O.riP.M.OT.2) 001 TO 50 
IF (R.oT.3.n) on in 20 
X=R*R/9, n 

IF (N.FO.l) on TO 10 
CALL AnAM{A,X,RJ) 

CALL AnAH(B,X,Y) 

BY=0. 5355l98.-:-Al.nr,(0,5=!=R )X:RJ+Y 
HZFRn=CMPLX( R.I,RY) 

IF (M.EO.O). return 
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10 CALL AnAM(C,XTY) 

BJ=R=i^Y 

CALL AnAM(n,X,Y) 

BY=n. <S366igR*ALnC{ 0.5*R )=I=RJ+Y/R 
HnWF=CMPLX( RJ,RY) 

RFTIIRM 
20 X=3.0/P 

IF (N.Fo.ii r,n in 30 
CALL AOAM(F,X,Y) 

FnnL=Y/soRT(R ) ; 

CALL AnAM(FtXtY) 

T=R+Y 

Rj=Fnni.*cns(T) 

RY=FOnL*SIW( T ) 

HZERn=CMPLX( RJ.RY) 

IF (N.FO.O) RETURN 
30 CALL AnAM(R,X,Y) 

FnnL=Y/SORT( R ) 

CALL AnAM(HtXtY) 

T=R+Y 

RJ=FnnL*COS(T) 

8Y=FnnL*SIM(T) 

HnME=CMPLXC RJ»RY) 

RETURN 

50 WRITE (-5,90) M.R 

90 FnRMAT( 32 H 0 SI CK DATA IN KAMKZl =J=OUIT=X N = , 1 2 , 2 X t 2 HR= , FI i « 3 ) 
CALL SYSTEM 
FNn 

SIIRROUTINE AnAM(C,X,Y} 

DIMENSION C(7) 

Y=X*C( 7) 
on 10 1=1,5 
10 Y=X=?-( C(7-I )+Y) 

Y=Y+C( 1 ) 

RETURN 

FND 

SURROUTINE FLIP( A,N,MI ,l,m,x,y,i at ) 

COMPLEX A(MI, 1) ,XM ) . Y{ 1 ) , D, R I OA , HOLD 
niMEMSION L(1),M(1) 

IF { lAT. GT.l ) 00 TO 150 

n=CMPLXC 1,0,0.0) 

no RO K=1,N 

L(K)=K 

M { K ) =K 

RIOA=A(K,K) 

no 20 J=K,N 

nn 20 i=K,M 

10 IF (CABS(3IGA).GE,CARS( A( t.J) ) ) RO TO 20 
Bir 7 A = A{I,J) 

UK) = I 
M(K)=J 




20 CONTIMUe 
J=UK ) 

IF (J.LE.K) r,n in 35 
nn 30 1=1 rW 

HnLO=-A( K t I ) 

A(K»I )=A( Jf T > 

30 A(j,i)=HnLn 
35 I=M(K) 

IF (I.LF.K) on TO 45 
nn 40 j=ifW 
HnLn=-A( .j,K) 

A( J»K )=A{ J, I ) 

40 A(J,I)=HnLO 

45 IF (CABS( BIOA) .NF.0,0) OR in 50 
n=CMPLX( O.OtO.O) 

RFTURM 

50 nn 55 1=1, M 

IF (I .FO.K) on in 55 
A{ I,K)=-A{ I,K ) /RIGA 
55 CnNTIMIIF 

nn 65 i=i,M 
nn 65 .)=1,M 

IF ( I.FO.K.nR, J.FO.K) on in 65 

A(I,J)=A( 1,K)^A(K,J)+A( ItJ) 

65 CnMTIMUF 

nn 75 J=1,M 
IF (j.fo.k) on in 75 

A(K,J )=A{K,J)/RIOA 
75 COMTIWIIF 

n=n*R I OA 

RO A(K,K)=1.0/RI0A 
K=M 

100 K=K-1 

IF (K.LF.O) on in 150 
I=L(K) 

IF { I.LF.K) on Tn 120 
nn no j=i,m 
H nin=A{ j,K ) 

A( J,K)=-A( J, T ) 

110 A(j,i)=:HnLn 
120 J=M(K) 

IF (J.LF.K) on Tn 100 
nn 130 T=i,M 
HnLn=A(K,i) 

A{K, I )=-A{ J, I ) 

130 A(j,i)=HnLn 
on Tn 100 
150 nn 200 1 = 1 , M 

Y( I )=CMPLX(0.0,n.0) 

nn 200 J=1,M 

200 Y( I )=A{ I, J)+Y( I ) 

RFTilRM 

FNn 

BLndK RATA 

CnMMnN/PIES/Pl,TPI,PlT,PIPTTYZTRFn,nTO 
nATA PI ,TPI , PIT, PIP I tYZ ,RF n,ni 0/3. 1415027,6.2031853, 
£1. 57079 63, 9. R 69 6044,0.0026525824,0.01745329,57.29578/ 
ENn VQ 


•C • INPUT FHRMAT fur PRHGKAM RICH 


SFPT, 1975 


* 4A> kC 4JU ^ 
» »«» **• »;• *»• *n I ^ 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C’ 

c 

c 

c 


CARn 1 FORMAT (1RA4) 


TITLF CARD: DSF UP TO 72 COLUMNS 


CARO 2 


FORMAT ( 12, I3t 5Fin,5) MORF,KOnE, ZFAC,WAVF, FIRST, LAST, INK 
MORF=0 THIS WILL RF THF LAST RUN FOR THIS OATA SFT 

MORE=l THERF are MORE DATA TO RF RFAO AFTFR THIS SET 

KOOF=0 COMPOTES RISTATIC SCATTFRIMO PATTERN 

KOOF=l COMPOTES R A CK SC ATTFR I MO PATTERN 

ZFAC A COMPLEX FACTOR MOLTIPLYINO ALL ELEMENT 

WAVE WAVFLFNGTH 

first INITIAL SCATTERING AMO INCIOENCE ANGLE 

last final amglf 

INK. ANGULAR INCRFMEMT 


CARO 3 


format ( I2,5X,7Fin.5) h , x A , Y A , XR , YR , ANG 
M 
I 

XA, YA, XR,YR 
ANG 


NOMRER OF SAMPLING POINTS ON THIS SEGMENT 
NORMALIZFO TNpFPANCE OF THF SFGMEWT 
SFGMENT ENOPniMTS 
angle SOBTFMOFn RY THF SEOMFNT 
ZFRO IN COL ? SHOTS OFF 

rearing of segments 


CARO 4 format ( 12, 13, Flo. 5) MORE ,K OOE , Z F AC , F I R ST , L AST , I MK 


THIS CARO IS OSEO ONLY 
ON CARO 2, MORF=l 


IF 


10 SPECIFICATIONS: 

?i:iMPOT OATA: f,:OOTPOT ( PR INTER ) : 7 : OUTPUT ( GR APH ICS ) 


C 
C 
C 
C 

C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

5 iJS ^ C 

c 
c 
c 
c 
=c 


C0MPLFX=:^R A(ino,lni),pHmQO),PINK{100),Z5(lO) 

C0MPLFX*R OFL. SOM, ZFAC 
RFA |.*4 1. AST. ink 

RFAL*A X{20n ) , Y( 200) ,XN‘( 2 00) ,YM( 200) ,S( 2 00) ,nSQ( 10) . ASUM{ 3A1 ) ,RSOM{ 3^1 ) 
INTFGFR*^ I0( 1«), LUMP (200) , I POL ( 2 ) , LV ( 1 00 ) , MM ( lOO ) 

DATA MH/lOn/ 

common /PIFS/ PI.TPI. pit, PIPI, Y7,RFn,niG 
OATA TPOL /'FE'=F«, » HHHH • / 

OATA IPP/2/ 

c RFAO TNPOT OATA AMO GFMERATF ROOY PROFILE 

5 RFAO (5,100) m 

RFAO (5,200) MORF, KOOF, ZFAC, WAWF, FIRST, LAST, ink 
WRITE ( 7, 100) in 
WRITE (7, 101) HAVE 
WRTTF(7, 19Q) FIRST, LAST, INK, IPP 
101 FORMAT!' L AM RO A= « , F5 .2 , ' MICRONS') 

1 Q 9 FORMAT! lX,3Fin.5, 13) 

IWHICH=1 

IF{REAL( ZFAC) .FO.O. ANO. AIMAG( ZFAC l.EO.O, ) 

SZFAC=(l.F-20,l.F-?0) 
write (6,150) in 
WRITE (6,300) 

CALL GEOf-U LUMP,X, Y. XM , YM, S » OSO, ZS , M ) 

MT=M/2 
LL=LUMP(M) 


original page JS 
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20 TF (KnnF.NF.O) P-n rn 
MIMC=1 

MRIT=1 + IFIX{ CLa.ST-FIRST)/IMK ) 
r,n TO 3n 
?H HRIT=n 

MINC = 1+IFI!<( ( LAST-FIRST) /I MK) 

30 WRITE ( A, 150) 10 

WRITE (ft, 400) IPOLI IPP) ,LL,MT,NIMCtMBIT,KAVE,ZFAC 
WRITE (ft, 42 5) FIRST 
3fl on 35 I=1 tLL 

ZS( I ) =ZS( I )=:=ZFAC 
35 nso( I )=nso( I ) /WAVE 
- XK=TPI/WAVF - 

C cnwSTRnr.T matrix flfmfmts 

IFdWWICH.Fn.?) on TO 37 

CALI. MTX(H,MH,XK,X,Y, XM,YM,nS0TLMMP,7S, A) 

C COMPlITF INCIOEMT FIELO AMO IN'VERT MATRIX 

37 TFTA=RFn*FIRST 
CT=CnS ( TFTA) 

ST=SIM( TFTA) 
no ftO I=2,M,2 

HOLOrr-XK* ( r.T4X( I )+ST4Y{ I ) ) 
nFL=CMPI_X( cns( HflLO) ,SIM{HDLn) ) 
ftO PIMK( 1/2 )=OFL*( XM( I )’-=CT + YM( n=^ST) 

CALL FLIP ( A,MT ,MH,L'/,mm ,P IMK ,PHI , I WH ICH ) 

C PRINT nilT CIIRRFMTS AMO FLFMFMT PRHPFRTIFS FOR FIRST AMOLF 

WRITE (ft, 500) 

IT=n 

on ft5 I=?,M,2 
IT=IT+I 

AMP=CARS ( PHI ( IT ) ) 

PHASE = niFi;4ATAM? (AIMAR(?Hl( IT)),REAL(PHI(IT))) 

ISFF=LIIMP( I) 

*65 WRITE (ft,250) I T, I SFR, X ( I ) , Y ( I ) , S ( T ) , OSO ( ISER ) , ZS { I SF(? ) , AMP » phASF 
c nnPF niiT thf appropriatf fiflo factors 

THF=FIRST-IHi< 

K=0 

IF (KOnE.eo.l ) GO TO 70 
WRITE (ft, 800) FIRST 
on TO 75 

70 WR I TE (ft, 600 ) 

75 THF = THF-^IMK 

IF (THF. GT. LAST) f;n TO 105 
IF { THE, EO. FIRST) GO TO 85 

C IM THF FOLLnwiMG LOOP, PINK IS MOT MECESSARILY THE INCIDEMT FIFLO 

TETA=RFn4THF 

CT=CnS(TFTA) 

ST=SIM( TFTA) 
on 80 J=2,M,2 

HnLn=-XK’:«(CT*X{ J) + ST*YIJ) ) 
nFL=CMPLX(cns(HnLO) ,siw{Hni.n) } 


75 


flO PrMK|J/?)=nFL-( XMC J)*CTh-YM( J)=!=ST) 

IF (KnnF.EO.O) r,n in 

CALL FLIP( LV,MM,priviK,PHI» ? ) 

R5 SUM=CMPLX( O^Ot n.o ) 

C ADD UP THE CURRFMTS 

nn 95 
JT=J/? 

9*5 SUM=SIIM+PHI( JT)*PIMK{ JT)’l'-nS0CLl)MP( J) ) 

SUMR=RFAL(SUM) 

.SUMT = AIMAfi( SUM ) 

S U M S 0 = S 1 1 M R =!: S U M R + S 1 1 M I * S U M I 
SCAT-10.^Al.nRin( 5I1MS0)+1.9A12 
K=K+1 

ASUM{K)=SUMSn 
RSUM{K )=REAL( SUM ) 

PHASE = ni(;*ATAM?{ SIIMI »SI)MR i 
WRITE! A, 901) THF, SCAT, PHASF 
WRITE(7»900) THF,SCATf PHASE 
on Tn 75 

10 5 n I FF= LAST-FIRST 

IFiniFF ,MF. 1R0,0 .AMO., RIFF .ME. 3A0.0) GH TH 205 
innn=Mnn(K,2) 

iFmiFF .en. BAo.n .amp, innn .me. d gh Tn 205 
IFCniFF .EO. IRO.O) FAC:^2.n 
IFiniFF .EO. SAO.O) FAC=1.0 
C*** SIMPSnw INTFGRATinM OF CRnSSECTinMS 
KLAST=K 

iFdnnn .ED. 0) klast=klast-3 

STGT=ASUM( 1 )- ^SUM( kLAST) 

STG=n.o 

nn 203 I=2»KL AST»2 
203 STG=S IG+ASliM ( I ) 

SIGT=SIGT4-4.0*SIG 

SIG=0.0 

KLAST=K-2 

nn 2in I=3tKLASTT2 
210 SIG=STG+ASUM( I ) 

STGT=SIGT+2.n^-RIG 

SlGT=SIGT=:aMK/3,0 

IFdnnn .ED. 0) SIGT=SIGT+3./R.^;IMK* 

K( ASUM( K-3 )+ASItM( K )+3.* ( ASUM( K-2 )+ASUM{ K-1 ) ) ) 
SIGT=SIGT*FAC=!=RFn/A.. 

pMf) SlMpsnM imtfgr ATinw 

KF=K 

IFIFAC ,F0. 1) KF=(K+l)/2 
STGF=-RSUM{KF) 

SrGA=-SIGT-fSrGF 
IF(KnnF) 220,220.240 
220 SrGTDR=10.0=!'ALnGLa{SIGT) 

IFISIGA .LF. 0.0) Gn Tn 225 
S IGADR-IO. 0* A LnG 10{SIGA ) 

Gn Tn 230 



225 WRTVB(6,R25) SIGT, Slf^At SIRF 
WRITE! Tta25) SIGTtSlGA,SIGF 
f 7 n tn 205 

230 5IGTnB=10.0*A|.nG10( Sir,!} 

WRITE! A, 850} SIGT, SIGTOR, SIGA^SrOAnR 
WRITE(7 tR 51) SIGTrSIGTOR^SIGAtSIGADR 
GO TO 205 

240 STGTnR=10.0*ALnG10( SIGT} 

WRrTB(4,aT5) SIGT.SIGTOB 

R25 FORMAT! ///,5Xt4nM*Fnil|.5!: MFOATIWE ARSORPTIOM CRnSSBCTION, //, 

£15X,5HSIGT=,FR. 5,5X,5HSIGA=,FB. 5T5X,5HSIGF=rFa.5 ) 

R50 FORMAT! /// , 1 5X , 5HS I GT= , FB . 5 , F8 .2 , 3H OR , 5X^ 5HS I G A= , F8 . 5 t FH . ?, 

£3H DR) 

R51 FORMAT!///tl5X,*SIGMA{T}/LAMRnA=»tFa.5fFa.2T’ OR^tSXt 
S'SIGHA! A}/LAMB nA=»tFa.5tFR,2, ' DP') 

875 FORMAT! ///» 15XT7HSIGAV/F=,Fa.5TFR.2t3W DR}- 
205 IF IMORF.EO.O) GO TO 5 

no 103 1=1, LL 
ZS! I }=ZS ! I ) /7FflC 
103 nso! I) =nSO( I IX^WAVE 
DFL=ZFAC 
pwAVE=WAVE 

RFAD !5,200) MORE , KOO F, Z FAC, WAVE? FT RST, LAST t INK 

WRTTB!7,100) id 

WRITE! 7, 101) waVF 

WRITEI7T19Q) FIRST, LAST, IMK,IPP 

IWHICH=1 

TF(REAL!ZFAC).EO.0.AWn.AIMAG(ZFAC) .Fu.n. ) 

£ZFAC=!1,F~20, l.F-?n) 

IF!REAL!7FAC ) ,FO,RFAL(nEL) . AWn.AIMAGIZFAC i.FO.ATMAG! OEL) 
£,AMD.WAVF,F0,PWAVF) IWHICH=2 
GO TO 20 

100 FORMAT ! 1BA4; 

150 FORMAT ! IHl, 1BA4 V 
200 format ! 12, 13, FF10,.5 ) 

250 FORMAT (?I5.7F10.5,Fin,3 } 

300 format (lOHOSEG M! iM , 1 1 X , 2*-HFWnP0IMTS OF THF SFGMEMTtl'lX, 

£1BHSEGMEMT PARAMFTFRS/llH MOM C BL LS , 6 X , 2HX A . BX , 2H Y A , BX , 2H XR , 8 X , 
£2HVR, AX. 24HAMGLF RAOIUS L FMGTH , ^X> 14HRF-Z IM-Z/) 

400 FORMAT I //31X, 14HKFY PARAMFTFRS// 

£16XT2lHlMCTnFMT POLAR rZATlOM, 22 X,1A17 
£laX,23WM{jM8FR OF SEGMFMTS IISFO,!?!/ 

£1 AX, 33HTnTAE WOMRFR OFV CELLS.^^^™^^ BODY, 111/ / 

£laXr35HMllMBFR OF IMCinFMT FIFLD niRBCT TOMS, 19/ 

£l6X,29Hi\illMRER OF RISTATIC D I R FCT lOMS , 1 1 5/ 
£14XtlOHWAVFLFMGTH,F34.5/ 

£IAX,16HrMPFnANGF F ATOR , FI A; 5 , ' 1 7 , fio . 5 
425 FORMAT !////, 3.5X,lRHSIiRFACF FlFLn nATA//27X, 

£29HF0R IMCIDFMT FIFLD OIRFCTIOM=, F7.? ) 

500 FORMAT (llHO I SBO ,4X ,4HX ! I ) ,OX r4HY ( I},6XT4HSt I) , 5X ,AHDS(K I ) 
£4X, 6HRS! I ) ,4X, (SHXS! I ) ,'^XT6HM0n( J ) ,^X, 6HARG! J ) / ) 






60n FORMAT ( lHlT27X»2RKRACKSCATTFRlNf; CROSS SECTinM//23X » 
g3f,HTHFTA S lOMA/LAMROA, OR PHASEtOFOI 

ROO FORMAT ( lHlt23X, 33HRISTATIC SCATTFRIMG CROSS SFCTinM/23X, 
S29HFDR INCinFNT FIFLO DIRFCTinM=T F7-2//23 Xt 
S 3AHTHFTA SIGM A/LAMRPA t OR PHASFtOFG) 

gOO FORMAT ( 15 XtF13.2» F14.2* F16-1 1 
901 FORMAT ( 1 5X , FI 3 . 2 1 F 14 . 2 1 FI S . 1 ) 

EMO 

SlIRROIITIME MIX ( M ,MH , XK .X T YtXWt VM » nSO f LUMPt ZSt A 1 
C RICH VERSION, H-POLAR I Z AT TOM 

REAL=i^4 X( 1) ,Y{ 1) ,XM( 1) .YN( 1 ) ,nSO(l ) 

IMTEGFR=!=4 LliMPC 1 ) 

CnMPLFX=!=8 A{MH»1) tZS( 1) THZTHlATHiR,n!lM,HlRLTHZA,HZRTHZBL 

COMMON /PIFS/ PItTPIt PIT, PIPI» YZ,RFD»niG 

IH=0 

IL=0 

no 10 1=2, M, 2 

IH=IH+1 

XI=X{I) 

YI=Y(I ) 

XMI=XN(I) 

YMI=YM{I} 

JH=0 

no 25 J= 2 ,M ,2 

JH=JH+1 

ISFG=IJJMP{J) 

ns=nso( IS EG) 

IF{IH.EO. JH)GO TO 50 

JA=J-1 

JB=J+1 

IF{ J.FO.M) JR=1 
XJA=X(JA) 

YJA=Y(JA) 

XJR=X(JP) 

YJR-Y(JB) 

XMJA=XM( JA)"'- 
YMJA=YN{JA) 

•XMJR=XM( JB) 

YMJR=YM(JR) 

RX=XI-X( J) 

RY=YI-Y{J) 

R-S0RT(RX*RX+RY-RY) 

SnSP = XMI=!-XW( J)+YNI^YN( J) 

RK=R^XK' - : 

CALL HAMKZK RK,0, HZ, nilM) . 

HZ=H2*S:nS'P.V:'^ ■ ■ 

RX=XI-XJR 

RY=YI-Y.IR 

R=S0RT( RX>:'RX+RY*RY ) 

SDR= ( RX*YWI-R Y=!iXN r r/R 

SnSP=XWI*XN.JR+YMT^YNJR 

RRK=R*XK 


78 








CALL HAMKZK RRKt2fHZR»HlR ) 

H1R=H1R*SDR 

H7.R=Sr)SP=;-HZR 

IF(XjA,FQ.XJRLcAMn,YJA,FD,YJRL. AMD, I.EO. IL ) DO TD 15 

RX=XI-XJA 

RY=YI-Y.IA 

R=SORT( RX*RX+RY*RY ) 

SnR=(RX':=YMI-RY*XMr)/R 
SnSP = XMI*'.--XMJA+YMI*YMjA 
RdK=R*XK 

CALL HAMKZK RAK,2,HZA,H1A) 

H1A=H1A^:'5DR 

HZA=SnSP'-:=HZA 

cn in 20 
V5 HlA=HlRL 
HZA=HZRL 

20 A( IH, JH) = PIT=i=nS7A,^{HZA+4.^i:HZ+HZR)-0.25^(HlR-HlA) 

XJRL-XJR 

YJRL=YJR 

H1BL=H1R 

HZRL=HZR 

on in 25 

50 A( IHt JW)=ZS1 rSFO)+ 

GCMPLX(PIT*DS/2.tl./ns/PlPI+nS=i=( ALnO(n.S )+0.02R79R37)/2. ) 

?5 IL=T 
10 COMTIMIIE 
RETIIRW 

FMO ■ 

SDRRniiTlWE r,EnM( L(lMP,XtY.XM,YM,StnSO»ZSTM ) 

C THIS VERSmM READS AMD aFMETATFSSEGMEWTS IM CnUMTER-CLnCKWISE DIRECTIUN. 
C THE SIIRFACF MUST RE CLnSFO. 

G PHINTS ARE CFMERFTFn AT THF START AMD MIDPniMTS DF EA ,H CELL; 

C THE START POIMT OF THF FIRST CFLL EVEMrUALLY CniMCIDE- 
C WITH THE FMD PniMT DF THE LAST CELL. 

CnMPLFX*R ZSClTtZ 

RF ALE'A X ( n , Y ( 1 ) , XM { 1 ) , YM { 1 ) , S ( I ) , DSD { 1 ) 

IMTFGFR^A LUMPtl 1 

CnMMnM /PIES/ PI,TPI, PIT, PIPI, YZ,REn,niG 

I=Q. . .. 

L=0' ", ■ ■ ■ ■ ■ ■ 

C.,,..BEAniMPMT parameters AMD PREPARE TO GEMERATE SAMPLIMG PQIMTS 
10 READ ( 5,200 ) M, Z , X A , Y A , XR , YR , AMG 
M. I T- 1 V on Tn 


OF POOR QUALrrf 
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Lim= 2 *m 

TX=XA-XB 

TY=YA-YR 

n=SQR TC TX-T X4-Ty^=TY ) 

L=L + 1 

ZS(L)=Z 

IF (AMG.EQ.O.O) Gn TO 
T=0.5-RED=!:AMG 


20 


TRX=T'/< + TY=’rCnTAM( T) 

TRY=TY-TX-C nl AM ( T ) 

RAn=-0.5*n/SIM(T) 

- ARC= 2 . 0 *RAn*T 
alf=t/m 

nSO( L)=2 AH^ALF 

GO in 30 

20 RAn=9F9.999 

ARC=n 

nson.)=n/N 
... START r-FMERATTMO 
30 DSSssO. 

ns=nso(L)/2. 
nn 50 JJ=1»L IM 
1 = 1+1 
J=JJ-1 
s(is=nss 
nss=nss+ns 
IJIMP(I)-L 

IF {AWg.fo.o.o} on in 4-0 

ARG=J=^^ALF 
EIMO=S IM( ARGl 
cns 0 =cns f ARG) 

X ( I ) = X A-0 . 5- ( T R ( 1 . 0-nn SO ) -TR Ys= S T MO I 
Yd) =YA-*-n. 5* ( TRXi'sS l.'in^TR Y* { 7 „ 0-Gn?Q ) ) 

XM( I ) =+0.5* { TRX*CnSO+TRY*SIMO ) /RAH 

YM ( I )= 0- 5* { TRX*SIMQ— TR Y*CnSO } /RAiO 
on TH SO 

40 Xd)=XA-0.5*J*TX/M 

Y( I )=YA-0.5*J*TY/m 
XNd)=“TY/0 
YM( I ) = TX/n 
50 COMTIMUF 

100 WRI TF ( 6. 300 ) Lf M, X A , YA. XR . YB » AMGfRAO, ARC T Z 
on TO 10 
120 M=I 
U=L 

200 FnRMAT|I?T3Xt7Fin,5) 

300 FnRMAT{ I 3 1 I Fr^-X , 4F1 0. 5 ♦ F7. 2 1 A-FIO- 5 ) 

RFTDRM 

FMO 

SUBRnilTTME HAMKZl t R tM.HZFBnTHnMF ) 

C, . . , .HANKFL FUMCTI HMS ARF of FIRST K IWn— J+TY 

c..... M=n RFTliRMS HZFRn 

G..... M=l RETIJRMS HOWP 

C..... M=2 RFTHRMS HZERn amp KOMF 

c stiRRniiTiNE reohirfs R>n 

C.....51IRRniiTIME APAM MUST BE SUPPLIFP RY MSFR 


DI MENS ION A(7),R{7)?C{T)tn(7)TEt7)tF(7]tG(7)TH(7) 
complex HZFRO^HnNF 

DATA A,R,C»n» E, F,r,*H/l. n,-?. 2499997, -0,31 A3flnA. ^ 

CO ,0444479 T-0 ,00394^4,0.0002 1 ,0 ,3n74ftft9l ^ 0.603 39346,-0.7433 0334 

SO. 233001 17, -0,042 612 14, 0.00427916,-0.00024046,0.5, -0.562499 A3, 

S0.21093373,-0.03P542R9, 0.004433 19,-0 .00031761,-0.00001109, 

S-0. 636619B, 0, 2212091. ?. 1632709,-1. 3164R27, 0.31 2 3951 ,-0.0400976 

50 .0027373.0. 7973343 6,-0.00000077,-0.00532 74, -0.00009512, 
sol 00137 2 37, -0.0007 2305, 0.00014476,-0.78539316,-0.04166397, 

S-0. 000039 54, 0.002 62573,-0 .00034125 ,-0 .000 7 9333 , 0 . 00M3 35 3 , 

50. 7973345 6. 0. 0^00015 6, 0.01639667, 0.0001 ( 1 05,-0 . 002495 11, 
SO.00113 633,-0.1. 0020033, — 2 .3 5619449 ,0.12499612,0.0000565, 

S-0. on 637879, 0,00074 34 8,0. 0007 9324, -0,00 0291 66/ 

IF (R.LF.0.0) on in 50 
TF (M.LT.O.nR.N.FT.2) GP TP 50 
IF (R. 07.3,0) on TP 20 
X=R^R/9.n 

IF (M.FO.I) r,n TP in 
CALL AnAM{A,X,BJ) 

CALL AnAM{R,X,Y) 

RY = 0, 63 66193=:= ALPOI n.5*R )=!=P.J+Y 
HZFRP=CMPLX( R.URY) 

IF (N.FP.O) RETURM 
10 CALL APAM(C,X,Y) 

RJ=R*Y 

call AnAM(n,X,Y) 

BY = 0 . 63 66193== ALPC i 0. 5=KR )=?B J+Y/R 
HPMF=CMPLX( B.I.BY) 

RFTPRM 
20 X=3.0/R 

IF (M.FO.l) on TP 30 

call ADANUF.X.Y) 

FnPL=Y/SnRT ( R ) 

CALL APAM( F,X,Y) 

T=R+Y 

RJ=FPPL*CPS { T ) 

RY=FnnL*SIP(T ) 

HZFRP=CMPLX( BJ .BY) 

IF (N.FO.O) RETURN 
30 CALL AOAM(C,X,Y) 

FPnL=Y/SPRT{R ) 

CALL AnAM{H,X,Y) 

j— ^ 

RJ=FnPL*CPS ( T> 

RY=FnPL^SIN(T ) 

HnWF=CMPLX( RJ,RY) 

RFTIIRM 

50 WRITF I 6,90) M,R 

90 FORMaTC 32H0.SICK DATA IN HAMKZl =!=DllTT* N=, 12 *2X , 2HR - 1 Fl I * 3 ) 
CALL SYSTEM 
FMn 

SURRniiTIME AnAM(C,X,Y) 
niMENSlPM C{7) 

Y=X*':=C(7) 





10 Y=x=s( C(7-I ) + Y) 

Y=Y+C( 1 ) 

RFTIIRM 

FMn 

SllRROMTIMF F L I P ( A ,M , M K |, , M , X , Y , I AT ) 

CrtMPLFX A(MI, 1 ),;<( 1 ) , Y( 1 ) , HtR JOAtHOLn 
niMFNSinM L( 1) .M{ 1 ) 

IF (lAT.OT.l) r,n TO 150 

n=CMpj.x( i.n.n.o) 
nn 80 K=i,M 

UKj^R 

M(K)=K 

RIGA=A{K,K) 

nn 20 J=K^W 
nn 20 i=KtM 

10 IF ( CARS ( RIGA ) . OF. CARS { A ( I » J ) ) } on TO 20 
RIGA=A(.I, J ) 

L(K)=r 
M(K) = .I 

20 COMTIWliF 
J=LtK ) 

IF (J, LF.K) GO Tn 35 
nn 30 I=1,M 
Hni.n=-A{K, I ) 
a{k,t)=a( j,n 
30 A(j,n=HnLn 
35 I=MC<) 

IF tl.LF.K) on Tn iS 
nn ^0 J=i,M 

MnLn=*A( 

A( J,K )=A( J, I ) 

^0 A(j,i)=HnLn 

A-5 IF (CARSIRIGA) .NF.n.n) GO TO 5n 
n=CMPLXfo.n,n.o ) 

RFTURM 

50 nn 5^5 T = l♦^• 

IF (I.FO.K) GO Tn 55 
A( I.K)=“A( NK )/RIGA 
53 CnWTIMItF 
nn A5 l=i,M 

nn 85 J= 1 ^M ORIGINM. PAGE IB 

IF { I.FO.K.nq. j,Ro,K ) Gn rn 85 OP POOR QUAIira 

An»J)=A{ I,K)*A('<,.i)+AtI,.J} 

85 cnNTlMni= 

nn 75 J=ItM 
IF (J.FO.K) Gn Tn 75 

A{K,J)=A(Kt J)/PIGA , . 

75 CnMTIMDF 
n=n*B I GA 

80 A(K,K)=1,0/RIGA 

K=M 
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mn K=K-i 

IF (K.LF.O) r,n Tn i 5 o 
r-UK) 

IF { I ,LF.K ) r,n Tn 120 
nn iin ,uri,M 
HnLn=A{ .UK ) 

A( J,K)=-Af.Un 
uo A(j»i)=HnLn 
I ?0 J=tM(K) 

IF (J.LF.K) r,n Tn ion 

nn 130 1=1, M 
HnLO=A(K,n 
AIK, I)=-;\( J,I 5 
130 A{j,n=HnLo 
on TO 100 
150 on 200 1 = 1, M 

Y( n-CMPLX(o.o,n,oj 
nn 200 .1=1, M 

.200 Y{ n=A{ I, J)*X{ J) + Y( I) 

RFTIIRM 

FMO 

BLOCK RATA 

CnMMni'i /P I E.S/Pl,TPI,PlT,PTPi,Y7»RFn.nTr; 

data pi , TPT.pt T.P TP T ,YZ,RFR. 010/3. 1A15°?7,6.2S31B5 3, 
fil. 57070 /.3TO,B6O.^0^-i,0.0n2n525R2A, 0. 01 7<i-5329T R7, 295 7B/ 
FMO 
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